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Abstract 

This paper studies robust transmission schemes for multiple-input single-output (MISO) wiretap 
channels. Both the cases of direct transmission and cooperative jamming with a helper are investigated 
with imperfect channel state information (CSI) for the eavesdropper links. Robust transmit covariance 
matrices are obtained based on worst-case secrecy rate maximization, under both individual and global 
power constraints. For the case of an individual power constraint, we show that the non-convex maximin 
optimization problem can be transformed into a quasiconvex problem that can be efficiently solved with 
existing methods. For a global power constraint, the joint optimization of the transmit covariance matrices 
and power allocation between the source and the helper is studied via geometric programming. We also 
study the robust wiretap transmission problem for the case with a quality-of-service constraint at the 
legitimate receiver. Numerical results show the advantage of the proposed robust design. In particular, 
for the global power constraint scenario, although cooperative jamming is not necessary for optimal 
transmission with perfect eavesdropper's CSI, we show that robust jamming support can increase the 
worst-case secrecy rate and lower the signal to interference-plus-noise ratio at Eve in the presence of 
channel mismatches between the transmitters and the eavesdropper. 
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I. Introduction 

Secure transmission in wireless networks is required in many applications. Traditionally, security 
is considered as an issue addressed above the physical (PHY) layer, and conventional approaches for 
ensuring confidentially are usually based on cryptographic methods. However, the broadcast nature of 
wireless transmission and the dynamic topology of mobile networks may introduce significant challenges 
to secret key transmission and management |fl], 0. Therefore, there has recently been considerable 
interest from an information-theoretic perspective in the use of physical layer mechanisms to improve 
the security of wireless transmissions. 

The theoretical basis of this area was initiated by Wyner, who introduced and studied the wiretap 
channel where the eavesdropper's received signal is a degraded version of the legitimate receiver's 
signal 0. The secrecy capacity was defined as the upper bound of all achievable secrecy rates, which 
guarantee that private messages can be reliably transmitted to the receiver and kept perfectly secret 
from the eavesdropper. Csiszar and Korner studied a more general non-degraded wiretap channel and 
considered transmission of secret messages over broadcast channels Q. Recently, considerable research 
has investigated secrecy in wiretap channels with multiple antennas ll5l- lfT3l . The secrecy capacity of 
the multiple-input multiple-output (MIMO) wiretap channel has been fully characterized in [8 ] and IfTTI . 
In particular, for multiple-input single-output (MISO) wiretap cannels, the optimal transmit covariance 
matrix was found to be single-stream beamforming obtained via a closed-form solution lfl4l . [fl~5l. While 
research in this area usually assumes global channel state information (CSI) is available at the transmitter, 
some other work has considered the case where only partial information of the eavesdropper CSI (ECSI) 
is available. The optimal transmit covariance matrix that achieves the ergodic secrecy capacity for the 
MISO wiretap channel was studied in |[T6l , where only statistical ECSI is assumed to be available. In IfTTI 
and [18], the problem was investigated from the perspective of maximizing the worst-case secrecy rate, 
where IfTTI considered the case that the eavesdropper's channel matrix is trace-bounded but otherwise 
unknown, and [18] assumed that the ECSI channel mismatch is norm-bounded. 

With the additional degrees of freedom in multi-antenna or multi-node systems, many papers have 
considered improving the secrecy rate through the use of artificial interference ||5), |fT3l , ||T9l , [20|. 
In these schemes, in addition to the information signals, part of the transmit power is allocated to 
jamming signals that selectively degrade the eavesdropper's channel while maintaining little interference 
to legitimate users. Some recent work has also considered using friendly helpers to provide jamming 
signals to confuse the eavesdropper GTI - BSI . This approach is often referred to as cooperative jamming. 
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In m , a cooperative jamming scheme is proposed for improving the achievable secrecy sum-rate for 
general Gaussian multiple access and two-way relay wiretap channels. The optimal transmit weights for 
multiple single-antenna helpers were studied in [24], [25 ], where a global power constraint was imposed. 
A similar case with individual power constraints was studied in (T7\ . In [28], the optimal beamforming 
strategy for a cooperative jammer was studied for the MISO case under a zero-forcing constraint that nulls 
the interference at the legitimate receiver. However, most of the previous work on cooperative jamming 
assumes perfect global channel state information, including CSI for the eavesdropper. This motivates us 
to investigate the case when the transmitters (including helpers) have only imperfect ECSI. 

In this paper, we study robust transmit precoder design for MISO wiretap channels with and without 
a helper. We assume that perfect CSI for the links to the legitimate user is available at both transmitters, 
while for the eavesdropper links there exist channel mismatches that are norm-bounded by some known 
constants lfl8l . |29l , lf30l . Following (9j, ifTTl . |[T5l . Gaussian inputs are assumed in the paper. We focus 
on obtaining robust transmit covariance matrices for 1) direct transmission (DT) without a helper and 
2) cooperative jamming (CJ) schemes with a helper (friendly jammer), based on maximizing the worst- 
case secrecy rate. Note that our work is different from [18] in that we consider robust transmission for 
not only information signals but also jamming signals. We begin by studying the optimization problem 
under an individual power constraint, and then focus on the more complicated case with a global power 
constraint. For the individual power constraint case, we show that the non-convex maximin problem of 
maximizing the worst-case secrecy rate can be converted into a quasiconvex problem that can be efficiently 
solved via existing methods. For the global power constraint case, we obtain the jointly optimal transmit 
covariance matrices and power allocation between the source and the helper via geometric programming. 
In addition, following |[T9l , |20ll . 11311 . we also consider the scenario where there is a quality-of-service 
(QoS) constraint at the legitimate node. In this case, we show that the introduction of the QoS constraint 
simplifies the optimization of the covariance and the power allocation. 

The organization of the paper is as follows. Section |II]describes the system model considered throughout 
the paper. In Section [TTll robust design of the transmit covariance matrix is studied for the direct 
transmission case. The robust cooperative jamming scheme is then investigated in Section |IyJ where 
both individual and global power constraints are considered. Section [Vj studies the case where a QoS 
constraint is required at the legitimate destination. The performance of the proposed robust transmission 
approaches are studied using several simulation examples in Section [VTl and conclusions are drawn in 
Section 17111 

The following notation is used in the paper: denotes expectation, (-) H the Hermitian transpose, 



January 19, 2013 



DRAFT 



3 



|| • || the Euclidean norm, (-)t the pseudo-inverse, tr(-) is the trace operator, IZ(-) represents the range 
space of a matrix, and I is an identity matrix of appropriate dimension. 



II. System Model 



We consider a MISO communication system with a source node (Alice), a helper (Helper), a destination 
(Bob), and an eavesdropper (Eve). The number of antennas possessed by Alice and the Helper are denoted 
by N a and Nh, respectively, while both Bob and Eve are single-antenna nodes. In this model, Alice sends 
private messages to Bob in the presence of Eve, who is able to eavesdrop on the link between Alice and 
Bob. The Helper can choose to be silent or to transmit artificial interference signals to confuse Eve. Both 
cases will be considered in the paper, and we refer to the former case as direct transmission (DT) and 
the latter as cooperative jamming (CJ). We assume that Alice and the Helper have perfect CSI for their 
links to Bob, but they have only imperfect CSI for their channels to Eve. We will consider cases with 
either individual or global power constraints. 

A. Direct Transmission 

When there is no support from the Helper, the received signals at Bob and Eve are given by 



where x is the signal vector transmitted by Alice, the covariance matrix of x is denoted by Q x = Ejxx^}, 
tr(Qai) < Ps where P$ is the transmit power constraint on Alice, and {h^, h e } are the 1 x iV a channel 
vectors for Bob and Eve, respectively. The terms n b and n e represent naturally occurring noise at Bob 
and Eve, and we assume that n b and n e are zero-mean circular complex Gaussian with variance a\ and 
<7g. We will assume without loss of generality that a 2 = a 2 = a 2 . 

B. Cooperative Jamming 

For the case where the Helper joins the network by transmitting an i.i.d. Gaussian interference signal 
z, Bob and Eve then receive 



y b = h b x + n b 



(la) 



y e = h e x + n e 



(lb) 



y b = h 6 x + g b z + n b 



(2a) 



y e = h e x + g e z + n e 



(2b) 
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where we denote Q z = K{zz H } and tr(Q^) < Pj. The cooperative jamming optimization problems 
that we consider in the paper will be subject to either an individual power constraint tr(Q x ) < Ps, 
tr(Q 2 ) < Pj, or to a global power constraint tr(Q z ) + tr(Q 2 ) < P. 

C. Channel Mismatch 

For the channels between the transmitters and Eve, only estimates h e and g e are available at Alice 
and the Helper, respectively. We define the channel error vectors as 

e h = h e - h e (3a) 

e 5 = ge - ge, (3b) 

and we assume that the channel mismatches lie in the bounded sets Eh = {e^ : [ 1 1 2 < e|} and 
£g = i e 9 • \\ e g\\ 2 — e <j}> wnere €h and e 5 are known constants. All the optimization problems in the 
paper are based on the e* h £ Eh and e* € E g that give the worst performance. 

III. Robust Direct Transmission 

In this section, we consider the scenario where there is no jamming support from the Helper. According 
to the signal model (Q]) and ([3]), the secrecy rate is ||9] 

hfeQ^hf \ - log 2 (l + 

A power constraint is imposed such that Q x € Q x = {Q x : Q x >z 0,tr(Q x ) < Ps}- For the case where 
perfect ECSI is available, the optimal Q x has been found to be unit-rank and the corresponding beam- 
former is the generalized eigenvector of the matrix pencil (cr 2 I + Psh^hb, a 2 I + Psh^h e ) corresponding 
to the largest generalized eigenvalue lTT5l . 

We consider the case where Alice does not have perfect knowledge of the channel to Eve, but only 
an estimate h e . We focus on optimizing the worst-case performance, where we maximize the secrecy 
rate for the worst channel mismatch in the bounded set Eh- Therefore, the optimization problem © 
becomes 

a 2 + hfcQ^hf 

max mm = = . (5) 

Q,g£L e h e£ h a 2 + (h e + e h )Q x (h e + e h ) H 

The difficulty in solving this problem comes from the inner minimization over e^. As will be discussed 

later, the minimization is actually a non-convex problem. However, we will show that through a proper 

transformation, problem (|5) can be converted to a solvable quasiconvex optimization problem. 



R s = log 2 ( 1 + -^h b Q x h^ - log 2 1 + — (h e + e h )Q x (h? + e h ) ) . (4) 
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Proposition 1: Problem © is equivalent to the following problem 

cr 2 + /ie 2 + tr[(Q :E + *)hfh e 



mm 



s.t. 



y o 



a 2 + tr(Q x h^h b 
* Qx 

Qx A*IjV - Qx 

tr(Qx) < P 5 
Q x h 0,/i > 0. 



Proof: The maximin problem in ([5]) can be transformed to 



max 



cj 2 + h b Q x h 



H 



s.t. <T + (hg + e h )Q x (h e + e h ) H < v 
e h eh < el, 

where the constraints can also be expressed as 

- e h Q x e% - 2Re(h e Q :E e^) - h e Q x hf - a 2 + v > 

- Chef! + e 2 h >0. 



(6a) 

(6b) 

(6c) 
(6d) 

(7a) 

(7b) 
(7c) 

(8) 
(9) 



>- 0. 



Using the 5 -procedure ll32l . we know that there exists an G C^" satisfying both the above inequalities 
if and only if there exists a \i > such that 

/A Na - Q x -Qx^e 

-h e Qx -h e Q x hf - a 2 - fiel + v 
Then we can use the property of the generalized Schur complement 11331 and rewrite ([TOl as 

a 2 + ne 2 h + h e Q x hf + h e Q x (^I Na - Q^Qzhf < v. 

Therefore, the maximin problem in d5) becomes 

a 2 + h b Q x h{? 



(10) 



(11) 



max 



Q,eQ^>o a 2 + fie 2 + h e Q x hf + h e Q x (fiI Na - Q s )+Q*hf 
which is equivalent to 



max 



a 2 + h b Q x h 



H 



cr 2 + ^ + heQ^h^ + he^hf^ 
s.t. Q x (nI Na - Q X ) ] Q X < *. 



(12) 

(13) 
(14) 
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Next, we use the Schur complement to convert (fl4T l into a linear matrix inequality (LMI), and the 
maximization problem is then given by 

a 2 + fie 2 + tr[(Q x + *)hf h e ] 



mm 

QxeQ x ,M>0,* 



s.t. 



>=0, 



(15a) 
(15b) 



<x 2 +tr(Q :c hfh 6 ) 

^ Qx fJ^N a ~ Qx 

which completes the proof. ■ 
Problem (O consists of a linear fractional objective function, which is thus quasiconvex, with a set of 
LMI constraints. Therefore, we can solve this problem efficiently via the bisection method [32). We first 
form the epigraph problem for (|6) as 



min t 



>- 



(16a) 
(16b) 

(16c) 
(16d) 



s.t. a 2 + n<? h + tr[(Q x + *)hf h e ] < t [a 2 + tr(Q ;E hf h 6 )] 

* Qx 

Qx ^N a - Q 

tr(Q x )<Ps,QxhO^>0, 
and the corresponding feasibility problem is then given by 
Find Q x ,/i,* 
s.t. a 2 + fie 2 h + tr[(Q x + *)hf h e ] < t [a 2 + tr(Q x hf h fe 

Qx ^N a — Qx 

tr(Q x )<P s ,QxhO,fi>0. 

Then we can use the following bisection method by solving the convex feasibility problem at each step 
until the interval that contains the optimal value converges: 



>- 



(17a) 
(17b) 

(17c) 
(17d) 



Algorithm III.l: Bisection method for obtaining Q x . 

initialize: set an interval [l,u] such that the optimal value t* € [l,u], and define a tolerance 5 > 0. 
while u — l> 5 

let t = (l + u)/2. 
do < Solve the feasibility problem in (fTTT ). 

if (fTTT ) is feasible, then u = t; else I = t. 
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For problem ©, the initial lower bound I can be chosen as any infeasible value, for example I = 
o- 2 +f^|ht,|| 2 , which is a lower bound since tr(Q x h^h&) < | [Qaj [ I2 • l|hfc|| 2 < Ps||hfc|| 2 . The upper bound 
u can be any feasible value, for instance, we can set u = a ^"^^[^"h^h" j hg ^ > wnere ^0 = -Ps> 
Qo = ftl, and * = QoW - Qo) f Qo- 

Note that the solution for the optimal covariance Q* obtained from Proposition [T] is already based on 
a hidden worst-case channel mismatch e* h . Next, we will explicitly express e* h under the norm-bounded 
constraint, which will be useful for the joint optimization in Section ITV-Bi The problem is formulated as 



max (h e + e h )Q*.(h e + e h ) H (18a) 



s.t. ||e ft ||<e ft . (18b) 

This is a non-convex problem since we want to maximize a convex function. However, we can still obtain 
the global optimum by solving its dual problem, as explained in the following proposition. 

Proposition 2: The worst-case channel mismatch for problem ([18} is given by e* h = h e Q* (AI — Q* X )K 
where A is the solution of the following problem 

max 7 (19a) 

A>0, 7 



s.t. 



y 0. (19b) 



AI - Q* Q*hf 
h e Q* -h e Q*hf -\e\- 1 
Proof: See Appendix lAl ■ 
Note that (fT9l is a semidefinite program (SDP) and hence can be solved efficiently using, for example, 
the interior-point method |32|. 

Thus, for the MISO wiretap channel, the problem of finding the beamformer for the worst-case 
eavesdropper channel can be converted from a non-convex maximin problem into a quasiconvex problem 
that can be solved using the bisection method. Both the optimal covariance matrix and the corresponding 
worst-case channel mismatch are obtained. In the following section, the MISO wiretap channel with a 
cooperative jammer will be investigated. 

IV. Robust Cooperative Jamming 

We now consider the case when the Helper provides cooperative jamming to improve the secrecy rate. 

According to the signal model in © and ©, the secrecy rate is 

„ , I, h t Q„h" \ , (, (h t +e h )QAk+e h ) H \ 

R ' = '° g2 I 1 + B Q.rf TP)- '° g2 I' + fe + ^jQM + »,)» TP)' <20) 
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We will first consider the optimization problem under individual power constraints, i.e., Q x € Q x = 
{Qx ■■ Qx h 0,tr(Q a; ) < P s } and Q 2 e Q z = {Q z : Q 2 >z 0,tr(Q z ) < Pj}, and then we investigate a 
more complicated case where a global power constraint tr(Q.j.) + tr(Q 2 ) < P is imposed. 

When the ECSI is perfectly known, maximization of R$ over both Q x and Q 2 requires an iterative 
search ll34l . However, it has also been shown that performance close to the optimal solution for the MISO 
wiretap channel can be obtained by simply forcing the helper's signal to be orthogonal to Bob's channel. 
Therefore, to simplify the solution when the ECSI is imperfect, we use a zero-forcing (ZF) constraint 
on the jamming signal for the CJ problem, which is equivalent to requiring that Q z gjf = 0. With the 
ZF constraint, the maximization of R$ with respect to Q 2 does not depend on Q x , although the optimal 
Q x still depends on Q 2 . Thus, we will first optimize Q 2 and then the optimal Q x can be calculated. 



A. Individual Power Constraint 

For the case of perfect ECSI, the optimal Q 2 under the ZF constraint is given by 

max g e Q 2 gf (21a) 
s-t- g 6 Q z gf = 0. (21b) 

Lemma 1: The optimal covariance matrix Q 2 G Q z for problem (1211 is rank one. 

Proof: See Appendix IB! ■ 

According to Lemma [U the optimal ZF jamming signal for the perfect ECSI case is also single-stream 
beamforming, i.e. Q 2 = Pj\vw H , where w is the unit-normalized one-dimensional beamformer for the 
Helper. We can rewrite problem (|2T1 as 

max |g e w| 2 (22a) 

w 

s.t. g fe w = (22b) 
w H w = 1. (22c) 
The solution for problem (|22l is referred to as the null-steering beamformer and is given by 

w - iid^-p^n (23) 

where P^^ = g^(g6g^)~ 1 g;> is the orthogonal projection onto the subspace spanned by gjf . The 
optimal information covariance matrix Q x , similar to the perfect ECSI case discussed in Section JIIJ 
is rank one and the corresponding beamformer is the generalized eigenvector of the matrix pencil 
P^h^he) with the largest generalized eigenvalue, where cr 2 — a 2 + g e Qzg,?- 
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For the case of imperfect ECSI, we still solve for the jamming covariance first, and the optimization 
problem becomes 



max min (g e + e g )Q z (g e + e g 
S-t. gbQzgb = 0. 



(24a) 
(24b) 



Proposition 3: Problem (1241 ) is equivalent to the following problem 

max tr[(Q 2 - *)gf g e ] - pe* 



s.t. 



>- 



(25a) 
(25b) 



Q z /uIat^ + Q 
tr(Q«) < Pj (25c) 
Q z h 0, /x > (25d) 

SbQzSt = °- ( 25e ) 

Proof: See Appendix O ■ 
Problem (|25T ) is an SDP that consists of a linear objective function, together with a set of LMI 
constraints. Therefore, we can solve this problem efficiently and obtain the optimal solution Q*. Note 
that although e g does not explicitly appear in (|25T ). the optimal robust covariance Q* is already based 
on the hidden worst-case e* that can be expressed explicitly through the following problem 

min (g e + e 9 )Q*(g e + (26a) 



s.t. 



|e s || < e g . 



(26b) 



This problem is similar to (fT8l i with the difference that (1261 is a convex problem and thus strong 
duality holds for (|26l i and its dual. The worst-case channel mismatch is provided through the following 
proposition. 

Proposition 4: The worst channel mismatch for problem (1261 is given by e* = — g e Q*(AI + Q^) -1 , 
where A is the solution of the following SDP problem 



max 7 

A>0,7 



s.t. 



ai + q: 



g e Q* g e Q*gf 



>- 0. 



(27a) 
(27b) 
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Proof: See Appendix |Pl ■ 
With solutions for Q* and e*, we can follow (H)-© and formulate the optimization problem over Q x 

as 

. ° 2 + (ge + epQ;(g e + e* g ) H + hfeQX 

max mm S- 1 , (28) 

Q,eS,e, G £ h o" 2 + (h e + e/ l )Q ;c (h e + e^) 11 

which can be solved with the same procedure as in Section JII] 



B. Global Power Constraint 

The previous sections only consider the design of the information covariance Q x and the jamming 
covariance Q z under the assumption that their individual power constraints Ps and Pj are fixed. This 
section investigates the joint optimization over Q x , Q z and the power allocation between Alice and the 
Helper, under the constraint that tr(Q x ) + tr(Q 2 ) = p\ + p2 < P. As with the previous case, we will 
assume a zero-forcing constraint for the helper's jamming signal at Bob. We will first optimize the power 
allocation for a pair of given Q x and Q z , and then we will provide an iterative algorithm that updates 
both the beamformers and the power allocation. 

First, for given and Q z , let Q x = piQ x and Q z = P2Q Z where Q x and Q z are normalized such 
that tr(Qa;) = 1 and tr(Q^) = 1. Hence the maximization of the secrecy rate (1201 ) with respect to p\ and 
P2 is equivalent to 

PiP2Cic 3 + picia 2 + P2C30- 2 + a 4 

max (29a) 

Pi,P2>0 p\C2 + P2C3 + o" 

s.t. pi+p 2 <P (29b) 

where 

ci =h 6 Q a; hf (30) 

c 2 ={h e + e h )Q x {h e + e h ) H (31) 

c 3 ={Se + e g )Q z (g e + e g ) H . (32) 

Since (I29ab is a quadratic fractional function, the optimization problem is hard to solve directly. However, 
we can use the single condensation method to solve this non-convex problem via a series of geometric 
programming (GP) steps [36]. GP is a class of non-linear optimization problems that can be readily 
turned into convex optimization problems, and hence a global optimum can be efficiently computed lf37l . 
Before applying the single condensation method, we give the following lemma 11361 : 
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Lemma 2: Given a posynomial 

m m 

/(x) = J>(x) = Y^^l^t ■ (33) 

i=l i=l 

where the exponents j5ji are arbitrary real numbers and a are positive constants, the following inequality 
holds: 

/(x)>/(x) = n(^) , (34) 
i=l V ai ' 

where aj > and YaLi a i = !• The inequality becomes an equality when a, = yj^y, in which case 
the monomial /(xo) is the best local approximation of the posynomial /(xo) near xo- 
Proof: We can rewrite /(x) as 

i=l v ' 

m / / ^ \ Q i 



i=i 



where (|36l ) holds according to the arithmetic-geometric mean inequality. When a« = jn^r . we have 
that «j (i = 1, • • • , m) are all positive coefficients with X]I=i a « = 1» an d the proof of equality is 
straightforward by inserting aj back into /(xo). ■ 
Using Lemma |2j let 

f(Pi,P2) = PiP2Cic 3 +p 1 c 1 a 2 +P2C3O" 2 + <7 4 , (37) 
and rewrite the numerator of (129a) as 

f(pi,Pz) =f(.Pl,P2) 



Q!l / 2\ a 2/ ^ 2\ a 3/4\ Q 4 

' (38) 



(\Oi/ 2\ "2 / 2\ a 3/4\ 

P1P2C1C3 \ / V\C\o A \ I P2C3C \ ( 
ax J \ a 2 J V «3 / \«4/ 



where 



"i = 77 r (39) 

f{Pl,P2) 

a 2 = 77 r (40) 

f{Pl,P2) 

a 3 = 77 r (41) 

f{Pl,P2) 

«4 = 77 r, (42) 
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so the optimization problem of (|29l l becomes 

Pic 2 +P2C3 + a 1 
mm = (43 a) 

Pi,p 2 >o hkpuvi) 

s.t. Pl + p 2 < P. (43b) 

The optimization problem stated above can be readily converted into the standard form for geometric 
programming problems, and (|43a) - (|43b) are posynomials of the GP problem. Therefore, the global optimal 
solution can be efficiently obtained for problem (l43l . Next, we can use the single condensation method 
to solve problem (|29l ), and the main steps are outlined in Algorithm II V. 1 1 



Algorithm IV. 1: Single condensation method for global power allocation 

initialize: power allocation and p^. 
iteration: 

For iteration A;, evaluate posynomial f( P i l \p^ ^) according to (l37T i. 



do < 



Compute aq according to (|39ll- (l42T ) using p^ ,p% ^ ■ 
Condense posynomial / into monomial / according to (|38l l. 
Solve problem (l43l with a single GP. 



(k) (k) 

Apply the resulting p\ and p 2 into the first step and loop until convergence. 



This successive optimization method is based on GP problems that can be solved using interior-point 
methods with polynomial-time complexity, and it has been proven in ll36l that the solution obtained 
using successive approximations for the single condensation method will efficiently converge to a point 
satisfying the KKT conditions of the original problem. 

Now we can conduct the joint optimization that considers both the information/jamming covariances 
and the power allocation between them. The main steps are outlined as follows: 

1) Initilize p\ = p 2 = -f . 

2) Let P s = pi, Pj = pi and solve problem dHJ, dHJ to obtain Q*„ Q*, e£ and e*. 

3) Let Q x = ;^qh), Qz = ^j=fj and solve problem (|29]> with Algorithm flV.il 

4) Apply the resulting p\ and p 2 to step 2 and loop until convergence. 

Since in each iteration the objective functions are maximized and the secrecy rate is increased, and 
since there is an upper bound for the secrecy rate (the case without an eavesdropper where the jamming 
power is zero), the procedure will converge to an optimum. Our extensive numerical experiments, some 
results of which are shown in Section [VT] further illustrate that the global optimum is obtained through 
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this procedure. Note that the non-robust counterpart of this procedure (where the CSI is assumed to be 
perfect) is similar but with the difference that instead of obtaining the covariance matrices from the robust 
maximin problems (1241 and (|28l l in step 2, Q* and Q* are the solutions of the non-robust problems in 

V. Robust Transmit Design with QoS Constraint 

In this section, we consider a slightly different problem in which there is a fixed constraint on the signal 
to interference-plus-noise ratio (SINR) of the legitimate link, so that the problem reduces to minimizing 
the SINR at the eavesdropper for the worst channel mismatch, under a given power constraint. As we 
will see, the addition of the SINR constraint actually simplifies the robust solution. 

A. Robust Direct Transmission 

When the Helper is silent, the optimization problem for the perfect-ECSI case is 

min h e Q x h^ (44a) 

Q,eQ, 

hf,Q x .hf 

s.t. =-2- > 7i , (44b) 

o~ 

where j t is the desired target SINR at Bob. When the channel mismatch is considered, the optimization 
problem is formulated as 

min max (h e + e h )Q x (h e + e h ) H (45a) 

Qx&Qx e h e£ h 

s.t. , b >lt- (45b) 

Comparing problem (l45l with ©, we see that (l43T l simplifies the fractional expression in © by intro- 
ducing an extra affine inequality constraint. Therefore, the procedure for solving © can be applied to 
solve (I43T ). and we have the similar result that (l45l is equivalent to the following optimization problem: 

min fie 2 h + tt[(Q x + *)hf h e ] (46a) 



s.t. 



h (46b) 



tr(Qx) < P (46c) 
Qx>=0,/i>0 (46d) 
trCQa-hf h 6 ) > a\. (46e) 
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Unlike the quasiconvex problem in Proposition [T] that requires a bisection method, (l46l ) is an SDP with 
a linear objective function, and thus can be solved efficiently. 

Note that by exploiting a similar form in (fTTT >. problem (|45T l is also equivalent to 

min fxe 2 h + h^hf + h e Q x (fj,I Na - QJ^^hf (47a) 

Q«eQ,.,/i>||Q»||2 

s.t. 5-2- > 7t . (47b) 



Denoting Q x = tr ^ x r , and based on the observation that h e Q x h,f > 0, heQa^/iljv,, — Qs)^Qa;h^ > 0, a 
relaxed solution for Q x of (l47l ) can be obtained by steering the eigenvectors of Q x to maximize hf,Q x h^ 
subject to a zero-forcing constraint Q x h^ = 0, while allocating the eigenvalues of Q x to satisfy the 
SINR constraint. Thus the relaxed optimization becomes 

max h fo Q x hf (48a) 

s.t. h e Q x h/ = 0, (48b) 

and Q* = - H Q*. Note that (l48l is also a relaxed optimization problem for its non -robust counterpart 
in d44l . because when there is enough transmit power to meet the SINR demand, Alice will choose a 
beamformer that will null the signal at Eve while aligning it as close to as possible. Therefore, in 
this case, the non-robust and robust covariance matrices, and thus the corresponding worst-case SINRs, 
are expected to be very close to each other. This will be illustrated by the simulation results of Section 

ED 

B. Robust Cooperative Jamming 

Unlike the previous section, in the cooperative jamming case with a QoS constraint, we will not need 
to introduce a ZF constraint on the jamming signal from the helper to simplify the problem. Note also that 
we only consider the global power constraint scenario, since the generalization of the proposed solution 
to the case of individual power constraints is straightforward. 

For the joint optimization of Q z and Q 2 under the QoS constraint, the optimization problem is given 

by 

(he + e h )Q x (h e + e h ) H 
mm max — — — (49a) 

q,,q, e hl e s (g e + e g )Q z (g e + e g ) H + a 1 

s.t. tr(Q x ) + tr(Q 2 ) < P (49b) 
| eft | < eft, |e s | < e g (49c) 
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h b Q x hf? 



> It- 



(49d) 



SbQzSb + <T 2 

Unlike the global power allocation problem discussed in Section IIV-B1 the QoS constraint in problem 
d49b simplifies the fractional quadratic expression in (|29l l into a linear fractional form with respect to Q x 
and Q z . Thus we can directly obtain the solution via the following proposition. 
Proposition 5: Problem (|49l ) is equivalent to the following problem 

lie\ + tr[(Q x + *)hf h 



s.t. 



mm 

Q.,q.,i»,*,w a 2 - ve\ + tr[(Q 2 - $)gfg e 

Qx mIat - Qx 

* Q, 

Q z uI Nh + Q z 

tr(Q x .) + tr(Q 2 ) < P 

Qx h o,Qz h o,fi> 0,1/ > o 

tr(Q«hf ht) > 7t(tr(Q 2 g?gt) + ^ 2 



>- 



>- 



(50a) 
(50b) 

(50c) 

(50d) 
(50e) 
(50f) 



Proof: See Appendix IE] ■ 
Problem (TSOb is also an SDP with a quasiconvex objective function, a set of LMIs and affine inequalities. 
Thus it can be solved via the bisection method discussed in Section HTH Note that for some cases when 
the SINR requirement at Bob can not be met for a given power constraint, the optimization problem will 
not be feasible and the transmission is assumed to be in outage. For the non-robust counterpart of the 
above QoS-constrained cooperative jamming problem, Alice will still use the solution of problem (144k 
and the Helper will use the remainder of the power for ZF jamming (same as problem (|2TI)). since the 
Helper has no information about the channel mismatch. 



VI. Numerical Results 

In this section, we present some numerical examples on the secrecy rate performance of the robust 
transmission schemes studied in the paper. For all examples, we assume Alice and the Helper both have 
four antennas, i.e. N a = = 4, while Bob and Eve each has one. The channel matrices are assumed 
to be composed of independent, zero-mean Gaussian random variables with unit variance. All results are 
calculated based on an average of 1000 independent trials. The background noise power is assumed to 
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be the same at Bob and Eve, a 2 = a 2 = 1, and the transmit power P is defined in dB relative to the 
noise power. 

We will examine the performance of the robust direct transmission (DT) scheme and the robust 
cooperative jamming (CJ) scheme under various power constraints, channel error bounds and QoS 
constraints. For purposes of comparison, we also examine the non-robust generalized eigenvector schemes 
(which will be referred to as GEV DT) discussed in Section HIT] , and the non-robust DT and CJ schemes 
discussed in Sections UTTTTVl 

Fig. Q] shows the worst-case secrecy rate as a function of transmit power under an individual power 
constraint, assuming P$ = Pj and the channel mismatch is e| = e 2 = 1.5. In general, when there exists 
channel mismatch, we see that the robust design for DT and CJ produces better performance in terms 
of the worst-case secrecy rate compared to their non-robust counterparts. Note that the non-robust CJ 
scheme, although possessing twice the available power as that for DT, performs even worse than the 
robust DT scheme for low transmit powers. This is due to the fact that the performance of the non-robust 
CJ scheme is degraded not only by the channel error between Alice and Eve, but also by that between 
the Helper and Eve. 

In Fig. [2l we compare the performance of the robust CJ scheme under both global and individual power 
constraints. In this case, we assume the global power limit P is lOdB, P$ and Pj are the individual power 
constraints for Alice and the Helper respectively, and P$ + Pj = P. The benefit of having the flexibility 
associated with a global power constraint over fixed individual power constraint is clearly evident. Also 
it can be seen that the proposed joint optimization procedure for the global power constraint achieves 
the optimal worst-case secrecy rate. When increases, a larger fraction of the transmit power must be 
devoted to jamming in order to reach the higher secrecy rate. 

The impact of the channel mismatch on the secrecy rate of the different schemes is presented in Fig. [3] 
The transmit power fraction for the robust CJ scheme is also plotted, and a global power constraint is 
used in this case. We assume P is 5dB, and the channel mismatch e| = e 2 . It can be observed that 
when the channel mismatch is zero, a jamming signal is not necessary, and all schemes achieve the same 
secrecy rate. However, when and e 2 increase, the robustness of the CJ scheme is more obvious, and 
the jamming fraction of the total transmit power also increases. Also it can be seen that the worst-case 
secrecy rate is much lower for the non-robust CJ scheme due to the fact that it is impacted by the channel 
mismatch from both the link between Alice and Eve, and the link between the Helper and Eve. 

Next, we consider examples for the case where a desired SINR constraint is imposed at Bob. P is 
assumed to be lOdB in these examples. In Fig. [4] the channel mismatch is given by e\ = e 2 = 0.5, 
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and we plot the measured SINR at Bob and Eve with an increasing SINR constraint at Bob. Since the 
SINR constraint is met in all cases, the curves for Bob all coincide. We see that in this QoS constraint 
scenario, the robust CJ scheme still shows a significant performance improvement by suppressing the 
SINR at Eve. An interesting observation is that the robust DT scheme for minimizing the worst-case 
SINR at Eve shows almost the same performance as the non-robust DT method, which is consistent 
with the analytical prediction discussed in Section IV-AI Therefore, the worst-case optimization for DT is 
unnecessary in this case. For the CJ schemes, when the SINR requirement at Bob is higher, more power 
will be allocated to Alice and less power will be used for jamming. Therefore, the performance of the 
CJ schemes approaches that of DT schemes. 

In Fig. [51 the impact of increasing e\ and e 2 g is depicted. In this example, the desired SINR constraint 
at Bob is set to be lOdB. We see that when no channel mismatch exists, all schemes give zero SINR 
at Eve while maintaining the QoS constraint at Bob, and in such cases the use of artificial noise (CJ) 
is unnecessary. Also we can observe that although the SINR at Eve increases for all of the methods as 
the channel mismatch becomes larger, the benefits of using robust cooperative jamming over the other 
algorithms is obvious. 

VII. Conclusions 

In this paper, we studied robust transmit designs for MISO wiretap channels with imperfect ECSI. 
Robust transmit covariance matrices were obtained for both direct transmission and cooperative jamming 
scenarios, based on worst-case secrecy rate maximization. For the case of individual power constraints, 
we transformed the non-convex optimization problem into a quasiconvex problem. For the global power 
constraint case, we proposed an algorithm for joint optimization of the transmit covariance matrices and 
power allocation. In addition, we also obtained robust transmit covariance matrices for the scenario where 
a QoS constraint is imposed at the legitimate receiver. The benefits of the robust designs were illustrated 
through numerical results. We conclude that although cooperative jamming is not helpful when perfect 
ECSI is available under a global power constraint, the worst-case secrecy rate can be increased and the 
SINR at Eve can be lowered by using jamming support from a helper when the ECSI is imperfect, 
provided that robust beamforming is employed. 
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Appendix A 
Proof of Proposition 2 



Problem (TT8T) can be rewritten as 



min - e h Q* x ef - 2Re(h e Q:,e£ ) - h e Q*h 

e h 

s.t. e h e%<e 2 h . 



(51a) 
(51b) 



This is a non-convex problem since its Hessian is negative semidefinite, i.e. — Q* < 0. The Lagrangian 
of this problem is 



L(e h , A) = - e h Q* x e% - 2Re(h e Q:ef ) - h^hf + A(e h e^ - <£) 
=e fc (AI - Q*)ef + 2Re(-h e Q*ef ) - h^hf - Xe 2 h 
where A > and the dual function is given by 

g(X) = mfL(e h , A) 



(52) 



e/, 



h e Q*hf - Xe 2 - h e Q* (AI - Q*)tQ*hf AI - Q* t 0, Q*hf e K(XI - Q*; 

otherwise 



(53) 



-oo 



where the unconstrained minimization of L(e^, A) with respect to is achieved when e/j = h e Q* (AI 
Q*)^. The dual problem is thus 



max - h e Q*hf - Xe 2 - h e Q:(AI - Q*)+Q*hf 
s.t. ai - q: y 0, Q*hf e K(XL - Q*). 
Using a Schur complement, the dual problem becomes the following SDP 



max 7 

A>0,7 



(54a) 
(54b) 



s.t. 



>- 0. 



(55a) 
(55b) 



AI - Q* Q*hf 
h e Q* -h e Q*.hf -Xe 2 "7 
Note that (I5T1) is usually called a frwsf region subproblem (TRS), and it has been proven that strong 
duality holds for TRS although the objective function is non-convex |[38l . Thus the optimal value of (|5"TT) 
and (1551) are the same. 
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Appendix B 
Proof of Lemma 1 

Since gbQzSb * s nonnegative, we rewrite the equality constraint as gbQzSb ^ £ < where e > is a 
arbitrarily small number. Denoting f{Q z ) = - q 1 h , the Lagrangian of problem (1211 is 

L(Q Z , A, 0, v) = f(Q z ) + A(tr(Q 2 ) - Pj) - tr(0Q 2 ) + vg b Q z gb (56) 

where A > 0, v > 0, and >z is the Lagrange multiplier associated with the inequality constraint 
Qz > 0. Note that the objective function in (|2TT ) is linear, and there exist strictly feasible points that 
satisfy both the inequality and equality constraints (e.g. Q 2 = J^-ww^, where w is a normalized vector 
orthogonal to g^). Therefore, according to Slater's theorem, the primal and dual optimal points of (1211) 
satisfy the Karush-Kuhn-Tucker (KKT) conditions 

tr(Q z ) < Pj, Q z > 0, gbQzgb < z (5V) 

A > 0,0^0 (58) 

tr(0Q 2 ) = (59) 

A(tr(Q a ) - Pj) = (60) 

e = -Td^? + XI + ^ U - <61) 

For the case that A = 0, according to (loTT i. we know that has a negative eigenvalue, which contradicts 
the fact that is positive semidefinite. Thus A can only be positive. For A > 0, according to d6TT > and 
since v > 0, has at least N — 1 positive eigenvalues, i.e. rank(0) > N — 1. 

Lemma 3 ( /l39l/): Given two N x N positive semidefinite matrices A and B with eigenvalues Aj(A) 

and Aj(B), respectively, arranged in non-increasing order, then 

JV 

tr(AB) > ^ A l (A)A 7V - 4 +i(B). (62) 

i=i 

Combining Lemma [3] and d59l , assuming Aj(0) and Aj(Q z ) are eigenvalues of and Q 2 , respectively, 
in non-increasing order, we have 

N 

tr(0Q z ) = 0>£ Xi(&)X N -i+i(Qz). (63) 
i=i 

Due to the fact that and Q z are both positive semidefinite matrices, we have 

N 

Y / \ t (e)\ N - i+l (Q z ) = 0. (64) 

1=1 
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Thus we can conclude that rank(@) ^ N, since otherwise all eigenvalues of Q 2 are zero and no jamming 
signals are transmitted. Combining this conclusion and the above observation that rank(@) > N — 1 , 
we can conclude that rank(@) = N — 1. Therefore, according to (164) , we have 

> 0, i = 1 
= 0, i = 2,--- ,N, 
which indicates that rank(Q 2 ) = 1, and the proof is complete. 



Aj(Q z ) { 



(65) 



Appendix C 
Proof of Proposition 3 

The proof is similar to the proof of Proposition Q] The maximin problem (1241 can be transformed to 

(66a) 



max v 



s.t. (g e + e 9 )Q 2 (g e + e g ) H > v 



The constraints in d66l ) can also be expressed as 

e g Q 2 ef + 2Re(g e Q 2 ef ) + g e Q 2 gf - v > 



e 9 ef + e 2 g > 0. 



(66b) 
(66c) 

(67a) 
(67b) 



According to the 5-procedure, there exists an e g £ satisfying both of the above inequalities if and 
only if there exists a fi > such that 



EeQz SeQzSe ~ ~ V 

Applying the Schur complement, we can rewrite (1681 as 



>- 0. 



tn S« 



Therefore, the maximin problem (I24> becomes 



max _ -/xe^ + g e Q 2 gf - g e Q 2 (^lAr h +Q 2 ) f Q 2 gf, 



QzeQ 2 ,/x>o 



which is equivalent to 



s.t. Q z (/jl Nh + Qz^Qz ^ *. 



(68) 

(69) 

(70) 

(71) 
(72) 
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Next, we use the Schur complement to turn the above constraint into an LMI. The maximization problem 
is then given by 

max - fiel + tr[(Q z - *)gf g e ] 



s.t. 



tr(Q 2 ) < Pj 



y o 



(73a) 
(73b) 

(73c) 
(73d) 



and the proof is complete. 



Appendix D 
Proof of Proposition 4 

The proof is along the same line as that in Appendix |A] Problem (l26l ) can be rewritten as 

min e ff Q*ef + 2Re(g e ef ) + g e Q*gf 



s.t. e g e g < e g , 



with the Lagrangian 



L(e g , A) = e 9 (AI + Q*)ef + 2Re(g e Q:ef ) + g e Q*gf - Xe 2 q 



(74a) 
(74b) 

(75) 



where A > 0. It can be verified that the minimum of L(e g ,X) with respect to e g is achieved when 



g e Q*(AI + Q*) , and hence the dual problem is 



max g e Q:gf - Xei - g e Q* z (XI + Q^^g 



s.t. AI + Q^O, Q*gf €K(XI + Q* Z ). 
Using a Schur complement, the dual problem becomes the following SDP 



max 7 

A>0,7 



s.t. 



ai+q: 



>=0, 



(76a) 
(76b) 

(77a) 
(77b) 



which completes the proof. 
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Appendix E 
Proof of Proposition 5 

The maximin problem (l49l can be transformed to 

v 

min - (78a) 

s.t. (h e + e h )Q x (h e + e h ) H < v (78b) 

v 2 + {ge + e g )Q z {g e + e g ) H >t (78c) 

tr(Q x )+tr(Q 2 ) <P (78d) 

\e h \ < e h (78e) 

|e 9 |<e g (78f) 

h*Q x hf = 7t (g 6 Q 2 gf + <r 2 ). (78g) 

Using the S-procedure similarly as in d8)-(H2) for (178b) . (|78e) and for (|78cT >. < T78fb . problem d78) can be 
rewritten as 

jue 2 + heQ^hf + h e Q x {^iI Na - Q^O^hf 
min — ( / y) 

Q-,Q.,/1>0,«/>0 CT 2 - l/e 2 + g e Q 2 gf - g e Q 2 (^IjV h + QzVQzge ' 

which is equivalent to 

/ie 2 + h e Q x hf + h e *hf 
mm — - (80a) 

s.t. Q x (^v a - Q^tQ, < * (80b) 
Q z (Mliv h + Q,) f Q, < * (80c) 
Using the Schur complement for the above constraints, the expression in (l50l can be obtained. 
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Fig. 3. Worst-case secrecy rate and transmit power fraction vs. channel mismatch, P = 5dB. 
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Fig. 4. Worst-case secrecy rate vs. SINR constraint at Bob, e| = — 0.5, P = lOdB. 
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